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Abstract 

We study free particles in a one-dimensional box with combinations of two types 
of boundary conditions: the Dirichlet condition and a one-parameter family of quasi- 
Neumann conditions at the two walls. We calculate energy spectra approximately and 
obtain equations of state having the same (one-dimensional) volume dependence as van 
der Waals equations of state. The dependence of the equations of state is examined for 
particles obeying Maxwell-Boltzmann, Bose-Einstein, or Fermi-Dirac statistics. Our 
results suggest that the deviation from ideal gas may also be realized as finite size 
effects due to the interaction between the particles and the walls. 
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1 Introduction 



In quantum mechanics, boundary conditions play an important role. A good example is 
Casimir force pQ predicted in 1948. It is also known that boundary conditions produce force 
in non-relativistic quantum mechanics [21 S] - In this paper, we focus on one-dimensional 
free particles and show that this force leads to equations of state whose one-dimensional 
volume or length dependence is the same as that of van der Waals equation of state. 

Consider a closed box. One may impose the conventional Dirichlet condition on wave 
functions at the walls of the box. In mathematics, we are, however, permitted to generalize 
the boundary condition in the context of self-adjoint extensions of Hamiltonian operators [5], 
where the Dirichlet and Neumann conditions are allowed just as special boundary conditions. 
Interestingly, it has been pointed out [6] that the generalized boundary conditions can be 
derived from the vanishing limits of widths of two step functions, which may be realized at, 
for example, junctions of two semiconductor layers. In this paper, we consider the Dirichlet 
condition and a class of boundary conditions including the Neumann condition, which we call 
quasi-Neumann conditions, to find how boundary conditions affect the physical property of 
the particles in the box. We show that such boundary conditions affect statistical quantities 
of finite size systems. Note that this effect has no relation with bulk properties of the system 
in the thermodynamic limit, where the effect that we will obtain disappears. 

Specifically, we approximately calculate the spectra of particles in the box by integral 
approximation and find that the spectra can lead to equations of state having the same 
dependence on length as van der Waals equation of state. Van der Waals equation was 
introduced by Johannes Diderik van der Waals in his Ph.D. thesis [7] and has been well- 
known to describe the behavior of real fluids. Van der Waals-like equations of state were 
derived from, for example, the Gaussian Markofnan process [8], the intermolecular interaction 
in classical statistical mechanics [21 EDI, the wall-molecular interaction in classical statistical 
mechanics [IT], or the electro- magnetic field [12J. Unlike van der Waals equation of state, 
our equations of state tend to that for ideal gas in the thermodynamic limit. In finite 
size system, the quantum boundary effects lead to similar equations of state as the above- 
mentioned effects. The free particles under the Dirichlet condition [13] and cyclic boundary 
condition [14] have been studied earlier, and here we consider combinations of the non-trivial 
boundary conditions and show that they admit equations of state whose length dependence 
is the same as that of van der Waals equation of state. 

This paper is organized as follows. In section 2, we briefly review the method of self- 
adjoint extensions to characterize the possible boundary conditions. In section 3, we study 
free particles in a one-dimensional box with the Dirichlet condition at its one wall and a 
quasi-Neumann condition at the other wall (see the left one of FigfTJ). We consider three 
types of statistics: Maxwell-Boltzmann statistics, Bose-Einstein statistics, and Fermi-Dirac 
statistics, for which equations of state are obtained. In section 4, we study free particles in 
a one-dimensional box with identical boundary conditions at its both walls (see the right 
one of FigflJ), where we assume Maxwell-Boltzmann statistics. In addition, we extend the 
one-dimensional box with identical quasi-Neumann conditions at both of the walls to a three- 
dimensional one. In section 5, we summarize our results and discuss the physical meanings 
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on the outcomes. 




Dirichlet q-Neumann q-Neumann q-Neumann 



Figure 1: The energy eigenfunctions and the energy levels in the boxes with different boundary conditions. 
The solid lines show the eigenfunctions. The height of dashed lines indicates the square root of the energy 
levels. 



2 On self-adjoint extensions of the Hamiltonian 

In this section, we give a short review on self-adjoint extensions based on |15j . Consider 
particles in a one-dimensional box with the length /. Each particle is governed by the 
Hamiltonian, 

h 2 d 2 

2m dx 2 ' ^ 

where x G [a, b] is a coordinate in the box. The Dirichlet conditions may be imposed 
on a wave function at the walls of the box. In general, infinitely large number of pairs 
of boundary conditions at both of the walls is permitted as long as the Hamiltonian is 
self-adjoint in mathematics. It is known that each pair of boundary conditions defines a 
distinctive wave function and spectrum. Thus, this generalization of boundary conditions 
has physical meaning. 

Let Dom(H) be the domain of the Hamiltonian. If if is a self-adjoint operator, then 
<pi,4>2 £ Dom(if) must satisfy 

b - h 2 d 2 , , , i r b . r h 2 d 2 



a 



2m I dx dx I 2m I dx dx 

Note that the above equations represent the conservation of the probability current density 
at both of the walls if d>\ = (bo. We assume that the wall of x = a is disconnected from 
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that of x = b, and for example, the cyclic boundary condition is excluded. The disconnected 
condition implies that each term of the second line of (12.21) must be equal to independently. 
We discuss only the boundary at x — a, since the boundary condition at x = b can be dealt 
with similarly. 

For convenience, we use <fit(x) introduced as 

^(*):=&(x)±iL^(s), (2.3) 

where L e M and i = 1,2. Note that L has the dimension of length. We can rewrite the 
condition (12. 2p as 

4>t*{a)<l>i{a) = fo*{a)fe{a). (2.4) 

The above equation means that the (one-dimensional) Hermite inner product of (f)f (a) and 
0^(a) is equal to that of (J)± (a) and </>j(a). Since fa and fa are arbitrary elements of Dom(/f), 
fal must be connected with fa by a (one-dimensional) unitary transformation specified by 
Dom(#): 

e ie <f>t(a) = 0r( a ). (2.5) 

Note that Dom(ff) is characterized by the unitary transformation vice versa. This procedure 
brings the condition of self-adjointness to a simple form: 

(e* e - l)faa) + iL{e %e + l)^(o) = 0. (2.6) 

ax 



This equation is known as Cheon-Fiilop-Tsutsui boundary equation[16j. 

We introduce Lg := — Lcot | instead of L and 8, then (12. 6ft can be written as follows: 

0(0 = 2^(0- (2-7) 

The above equation contains the arbitrary parameter Lg. In mathematics, the Hamiltonian 
( 12. ip is not a bounded operator; therefore it cannot be defined in the entire Hilbert space 
but a dense subspace of it. The domain of the Hamiltonian cannot be specified uniquely, 
either. The arbitrary parameter Lg implies this ununiqueness. In particular, Lg = and 
Lg = oo imply the Dirichlet condition and the Neumann condition, respectively. Note that 
the boundary condition (12.71) has been derived [5] from the limits of two step functions, 
whose ratio between the heights and widths characterizes Lg 



3 The box with the Dirichlet and a quasi- Neumann 
conditions 

In this section, we study the box with the Dirichlet and a quasi-Neumann conditions. 
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Fi gure 2: The thin lines represent y = cot kl and the dashed lines show their linear approximation. The 
thick line is y = j^Yr ^ ne cross i n g points of the thick and thin lines yield the solutions of (|3.2[) . Those of 
the thick and dashed lines yield the approximation solutions. 

Let x G [0, 1] be a coordinate in the box. We assume that our wave function satisfies 
the Dirichlet condition at x = and a boundary condition characterized by Lq at x = I. 
We restrict ourselves to \Lq\ > I at x = I, which is required technically to ensure the 
validity of our approximation. A class of such boundary conditions includes the Neumann 
condition, but not the Dirichlet condition; accordingly, we call such a boundary condition a 
quasi- Neumann condition. 

A wave function of each particle is given by 

<f>(x) = Asinkx, (k G R), (3.1) 
and its spectrum condition is obtained by 

-f- = kl cot kl. (3.2) 

Note that the above equation is a transcendental equation; therefore we cannot solve it 
explicitly. 



3.1 Calculation of approximated spectrum 



Now we study the equation (13. 2p and obtain its approximated solutions. The spectrum 
condition (13.21) has a solution for each region, (n — l)n < kl < wr, where n G N. In addition, 
the solution kl tends to 7r(n — |) in the limit n — > oo (see Figj2]). We expand cot kl as 



cut kl = ( n — - )rr — kl -\ O 



1 \3 

kl — nn H — tx 
2 



(3.3) 
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For each region, the following inequality is fulfilled: 





< cot kl = 


I 1 


(n "W — kl 




\ 2) 




La kl 



(3.4) 



This inequality leads to 



[n W — kl 

v 2 1 



< 



L e7 i(n-\) 



O 



(3.5) 



The approximation solution of (13. 2p for each region, {n — 1)tt < kl < rnr, is therefore obtained 
by 



kl = 7r( n 



I 



2/ -nL e {n-\) 



I 



Len(n- \) 



(3.6) 



where the error term is estimated from (13. 5p . The solutions lead to energy spectrum of each 
particle E n as^| 



E„ 



K 
P 



I \2 1 



1\2 V 



n 



1-0 



Lon(n- |) 



k{u-\) 2 v 



where 



h 2 



v :- 



2m ' itlLq 
and we ignore the terms of the order by our quasi-Neumann condition. 



(3.7) 



(3.8) 



3.2 Statistical quantity and equation of state 

In this section, we calculate statistical quantities and derive equations of state from the 
spectrum in high temperature region ^ < 1. We assume iV particles obeying Maxwell- 
Boltzmann statistics, Bose-Einstein statistics, and Fermi-Dirac statistics. 

2 The approximated spectrum (|3.7p and (|4. 14[) were studied in a more general context: e.g., [T7], where 
the approximation is valid in the region n > no for some no. Our formulae (13. 7p and (|4.14j) are valid for all 
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3.2.1 The case of Maxwell-Boltzmann statistics 



Now, we study the case of Maxwell-Boltzmann statistics. By using (1A.4[) . we obtain partition 



function Z c as follows: 



^exp(-/3£ n ) \ 

n=l ' 

= 6XP (^r){/ dxex P (-^ 2 )| 



A' 



-J— exp(-f-). (3.9) 



2V fJnJ ^\ I 



Let p c be force acting on the walls of the box. We can, based on statistical mechanics, derive 
p c from following equation: 

1 dZc to m\ 

Vc = Qr7 J7 ■ (3.10) 



(3Z C dl 

The equation of state in this case is thus given by 

Nu\ , N 



The given equation differs from ideal gas law in the term We refer to this term as a 
force correction term. Note that (13. lip becomes the equation of state for ideal gas if the 
wave function satisfies the Neumann condition at x — I. This implies that equation of state 
for ideal gas cannot be derived under a pair of Dirichlet conditions, which will be discussed 
in section 14. 11 

The equation of state (13.1 ip resembles van der Waals equation of state: 

Pvdw H p J [l- VvdwN) = —. (3.12) 

Van der Waals equation of state is different from ideal gas law in the terms N ^ dW and 
— Cvdw^V- We call the former a force correction term and the latter a length correction term. 
Van der Waals equation differs from (13. lip in two points. One is that van der Waals equation 
of state has the length correction term, while (13.111) does not. In section HI we will derive 
equations that possess length correction terms, which are however different from van der 
Waals one. The other is that the force correction term in (13. lip is in proportion to N, while 
that of van der Waals equation of state is in proportion to iV 2 ; therefore, the force correction 
term in (13.111) vanishes in the thermodynamic limit I — > oo under the constant density. This 
is reasonable because the boundary effects are expected to disappear in the thermodynamic 
limit, where the property of bulk gas becomes dominant. Our equation is thus meaningful 
only in finite length systems. We will further discuss this point in the last section. 
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3.2.2 The cases of Bose-Einstein statistics and Fermi-Dirac statistics 



Now, we study the case of quantum statistics. Let bosonic variables be subscripted by + 
and fermionic variables by — . 

We express, by using integral approximation (1A.4[) . the particle number N as a function 
of chemical potential [x± like that: 



1 



^ exp {(3(E n -/i±)} =F 1 



i nr 



±exp|/?(/i ± + j) 



(3.13) 



where Li s (x) is polylogarithm function (see Appendix [B]). Note that the integral approxima- 
tion of the bosonic case is violated in the region /?(//+ + f) ~ since the integrand becomes 
infinity at /?Ut+ + j) = 0. We will discuss this case in section 13731 

Analogously, the force in this case p± is obtained by 



DC 



dE n 



P± E di exp{f3(E n -^)} T l 

oo 



\± V 1 2) I 2 

< exp [/?(£„ -ft-)} T 1 



2pV k(3 a 



±exp<J/?(y + /i±) 



iVi/ 
I 5 "' 



(3.14) 



We can rewrite the equation (13. 13j) as 



P 



Lh{±e^ ±+ f)}j 2 . 



(3.15) 



Since Lii(±e y ) is a monotone function of y, there exists the inverse function; therefore, 

/3(/i+ + j) can be regarded as a function of K ^p . For this reasons, we introduce R± ^!^- 
by 



k(JN 2 

~p~ 



Li 3 

2 


± exp < 


[/5(f 


+ M±) 


n 


Lii 

2 


± exp < 


[/5(f 







(3.16) 



whose behavior is shown in Fig. [3J 

From the above equation, we can rewrite (I3.14p as 



Nu\, N /Kf3N 2 

p± + — )l = —R± 



P ) 



P 



(3.17) 




k/3N 2 



R 

6 

4 

2 
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k/3N 2 
I 2 



Figure 3: These figures show the behavior of R±. (Left) Bosons. It shows that R+ is a monotonically 
decreasing function of . (Right) Fermions. It shows that R- is a monotonically increasing function of 

0. 



The graphs imply that both R± tend to 1 in the limit 



We regard the above equation as the equation of state in this case, which also has the term, 
as is the case for (13. lip . We call this term a force correction term as well. In the 
bosonic case, a similar equation without the force correction term was derived from the 
cyclic boundary condition in [T3j. Note that the equation (13.111) is derived from replacing 
R± with 1. 



3.3 Behavior in the limit / — > 

In the previous section, we calculate the statistics quantities in the high temperature region 
w < 1. In this section, we study the behavior of the equations of state in the limit /3 — > oo. 



3.3.1 The case of Maxwell-Boltzmann statistics 



In the region j§- 3> 1, the integral approximation is not better than first term approximation: 



Then p c is given by 



Z c = exp(-^^ + — — ). (3.18) 



(Pc + -^ = ^. (3.19) 



We call the term ^ a force correction term as well. Note that the term is the same shape 
as those of fEUTj) and fl3TT?j) . 



3.3.2 The case of Bose-Einstein statistics 



In the bosonic case, the integral approximation is not verified in the region + y) ~ 0, 
since the integrand of (IA.4I) diverges at + j) =0; therefore the second term of the right 
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hand side of AA.4D cannot be neglected. In this region, we use the following approximation: 

1 



^ 1 e X p{%(n-iy-l3(j i+ + z)}-l 

1 f°° dx 



exp{§-/3(/i + + f)}-l Ji exp{^ar»-/?(Ai + + f)}-l' 
An upper bound of the integral term is given by 



(3.20) 



dx f°° dx 

< 



i exp{^x 2 -/^ + + f)}-l Ji $ x *-p(n + + !t)-l 

<£• (3 - 2i) 

If y 3> the first term of the second line of (I3.20p is much bigger than the integral 
term. In other words, the average number of the particles in the ground state is much bigger 
than that in the excited states. Bose-Einstein condensation therefore happens in that region. 
Note that the condition, y 3> 4^, is independent of the boundary parameter L e . In this 
region, we can neglect the integral term of (I3.20p ; therefore the equation of state is given by 
as follows: 

(P+ + -p) l =2p- ( 3 - 22 ) 
The above equation is identical to (13 . 19[) based on Maxwell-Boltzmann statistics. 



3.3.3 The case of Fermi-Dirac statistics 

In the fermionic case, the equation of state (13.171) is still valid in the limit of j3 — > oo. From 
appendix [B], the asymptotic behavior of R-(z) in the limit of z — > oo is 

itL(z) -^-z. (3.23) 

The equation of state thus tends to 

(P- + ~jr) 1 = -gp-- ( 3 - 24 ) 

We also call the term ^ a force correction term. Note that the force correction term is the 
same shape as those of fl3TTTD . fl3TT7j) . (|515|) . and <KW\i . 



4 Pairs of identical boundary conditions 

We have studied the statistical quantities of the system in the box with the Dirichlet condition 
and a quasi-Neumann condition in the previous section. In this section, we assume pairs 



10 



of identical boundary conditions: a pair of the Dirichlet conditions and a pair of identical 
quasi- Neumann conditions. In addition, we extend the one-dimensional box with a pair 
of quasi-Neumann conditions to a three-dimensional box. We consider particles obeying 
Maxwell-Boltzmann statistics since the quantum statistics cases are hard to study by using 
the method developed in section [3][1 



4.1 One-dimensional box with the Dirichlet condition 

In section 13.2.11 we have shown that the equation of state under Maxwell-Boltzmann statis- 
tics with the pair of the Dirichlet condition and the Neumann condition is identical to that 
for ideal gas. In this section, we show that the equation of state of the pair of the Dirichlet 
conditions is different from that for ideal gas. 

We impose the Dirichlet condition at the walls of the box. Let E n be energy of each 
particle given by 

2 



En = — , (4.1) 



K71 
1 

where n G N. By using flA.8j) . we obtain partition function Z DD as follows: 



z DD = {J2 e M-PE n )} 



N 



n=l 

dxexp ( - —x ) - -j 

/ nr 1 \ n 

-J— - -) . (4.2) 
.2\j k(3 2 J K ' 

Let p DD be the force acting on the wall. The equation of state in this case is thus given by 

We call the term —\J~^- a length correction term. The above equation of state is clearly 
different from ideal gas law, which is derived from classical statistical mechanics; therefore 
quantum effect is responsible for this difference. 

Tsutsui et al. have studied a one-dimensional box with a partition at the center of the 
box [21 E] - They imposed the Dirichlet conditions on wave functions at the walls of the box. 
They also imposed the Dirichlet condition at one side of the partition and the Neumann 
condition at the other side of the partition. It was shown that the partition is, in high 
temperature limit /? — > 0, subjected to force Ap given by 

n 



3 The particle number is not a monotone function like (|3.13[k therefore, we cannot repeat the discussion 
of section |31 
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where I is the length of each half of the box. 

Below, we derive Ap from the equations of state (13.111) and (14. 3p . In the high temperature 
limit, the length correction term is smaller than l; therefore we approximately rewrite p DD 
as follows: 

N 1 



P DD 




(4.5) 

where p DN is the force of Maxwell-Boltzmann particles in the box with a pair of the Dirichlet 
condition and the Neumann condition. Since the difference of the forces p DD — p DN is given 
by Ap, our results agree with that of [21 [3] in the high temperature limit. 



4.2 One-dimensional box with a quasi-Neumann Condition 

In this section, we study particles in the box with a quasi-Neumann condition. Let x 6 [— -, | ] 
be a coordinate in the box. We assume identical boundary conditions at the walls. In 
accordance with [18], the boundary conditions at x = | and characterized 
by the boundary parameters Lg and — Lg, respectively. The sign change of the boundary 
parameter is caused by reversing the direction where the particles collide on the wall. The 
wave functions are given by three the types: 

<p = cos kix (4.6) 
(j> = sin k 2 x (4-7) 
4> = cosh k 3 x, (4.8) 

where fci, k 2 , and k 3 satisfy following quantization conditions: 

I , , k\l , 

A;i/tan— (4.9) 



L e 2 

1 1,1 ^ 

— = kol cot — 

L e 2 

— = k 3 l tanh - 
U 2 



/ k 2 l cot M (4.10) 
k I 

/c 3 Z tanh— . (4.11) 



The approximation (13. 3p and 



tanx = x - mr + 0[(x - hit) 3 ] (4.12) 

tanh x = x + O [x 3 ] (4.13) 

lead to energy spectrum of each particle as 

n(n-l) 2 2v 
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where n G N. By using equation (IA.8I) . we obtain partition function Z NN as follows: 

oo 

z nn = {x>p(-^)r 



N 



n=l 



exp 



/ 2Nv(3 
\ I 




dxexp ( - ^V) + ^} 



/2Nup\/i nr l^ JV 

exp ( — ; — ) ( - J— + - 



V I J\2\np 2, 
Let force acting on the wall be p NN . An equation of state in this case is given by 



(4.15) 



p *» + ^( l+ IP* 




N 
J' 



(4.16) 



We call the term a force correction term and the term y ^ a length correction term. 

Note that v = implies the Neumann condition. Let p N be force under the pair of the 
Neumann conditions, the equation of state is given by 



v N (i+\r- 

V 7T 



N 

J' 



(4.17) 



The above equation is different from the equation of the pair of the Dirichlet conditions 

flop. 



4.3 Three-dimensional box with a quasi- Neumann condition 

Here, we study a three-dimensional box with a quasi- Neumann condition at the wall of box. 
Let (x, y, z) be a three-dimensional Cartesian coordinate. Consider a cuboid box whose 
median point is the origin of the coordinate. Let its sides be parallel to the coordinate axes. 
We introduce l x ,l y , and l z as the lengths of the sides. As is the previous section, we impose 
identical boundary conditions at the walls of the box: the boundary condition at x = y, 
y — ^, and z = l f are characterized by L e and the boundary condition atx = — y = — y, 
and z = — l -f are characterized by — Lg. We introduce quantum numbers for x, y, and z 
directions G N, which characterize energy of a particle, E nx)n nz} as 

((n x -l) 2 K-l) 2 (n z -l) 2 ) (\ 1 1\ 

*„„„„„, = «{^- + LyL + } " 2 "U + T, + r)' < 4 - 18 ' 

The energy implies that force is different by the directions. We introduce p x as force of the 
x direction, the equation of state is given by 

+ V + 4 V7) = ? (419) 
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where S x = l y l z and V = l x S x . The force of the other directions is obtained by the same 
calculation. Note that the terms showing the aeolotropic behavior of pressure vanish in the 
thermodynamic limit V — > oo under the constant density condition, y = constant. The 
equation of state (I4.19p thus predicts no aeolotropic bulk gas but shows the dependence of 
boundary effects on the scale of the system. 

5 Conclusion and Discussion 

We studied free particles in a closed box whose boundary conditions are generalized by 
self-adjoint extensions of the Hamiltonian. We solved the energy spectrum conditions ap- 
proximately and obtained the equations of state from the energy spectra by integral ap- 
proximation. Note that our results are significant primarily for finite size systems since the 
boundary effects disappear in the thermodynamic limit; therefore, our equations of state do 
not describe the bulk properties of free particles at the thermodynamic limit. 

We assume three different combinations of the boundary conditions. First, we consid- 
ered the one-dimensional box with the Dirichlet condition and a quasi-Neumann condition. 
We assumed three types of statistics: Maxwell-Boltzmann, Bose-Einstein, and Fermi-Dirac 
statistics, and obtained the equations of state (13.111) and (13.171) and evaluated the behavior 
of the equations in the limit I — > 0, arriving at (13.191) . (13.221) . and (13 .24p . Although these 
five equations are different in statistics and in the region of /, they share the same force cor- 
rection terms. The common force correction terms have the same dependence on length as 
that of van der Waals equation. The force correction terms are dependent on the boundary 
parameter, whereas other terms are independent. 

Second, we studied Maxwell-Boltzmann particles under the Dirichlet condition at both 
of the walls. We found that the equation of state (14. 3 p is no longer that of ideal gas due to 
the length correction term in (14.31) . which is different from that of van der Waals equation 
of state (I3.12p . It has been studied in [HE] that there is a difference between two pressures, 
one of which obtained under the pair of the Dirichlet conditions and the other evaluated 
under the pair of the Dirichlet and the Neumann conditions. The equations of state of 
(13. lip and (14 .3p lead to the high temperature limit of the pressure difference (14.41) . which is 
an increasing function of temperature and vanishes in the thermodynamic limit. 

Third, we studied particles obeying Maxwell-Boltzmann statistics under a quasi Neumann 
condition at both of the walls. We obtained the equation of state (14.161) . which has both the 
length correction term and the force correction term. The strength of the force correction 
term in this case is twice of that of the pair of the Dirichlet condition and a quasi Neumann 
condition. The length correction term in this case is similar to that of the pair of the Dirichlet 
conditions except for its sign. Note that the equation of the pair of the Neumann conditions 
(I4.17P is different from that of the pair of the Dirichlet conditions. 

These three results clearly show that boundary conditions can yield force correction terms 
that have inverse square dependence on the length, analogously to van der Waals equation. 
However, the force correction term in van der Waals equation and those of our equations are 
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different in particle number dependence. Our equations of state are influenced by boundary 
effects which vanish in the thermodynamic limit, and hence they are meaningful only in the 
finite length system. In the case of van der Waals equation, the order iV 2 is regarded as the 
number of the inter-molecular interactions, which is N(N — 1) ~ TV 2 since the iV particles 
interact with each other. On the other hand, the force correction term of the pair of the 
Dirichlet condition and a quasi Neumann condition is in proportion to iV and that of the 
pair of quasi Neumann conditions is twice of that of the pair of the Dirichlet condition and 
a quasi Neumann condition. There is no force correction term in the equation of state of 
the pair of the Dirichlet conditions. Interestingly, the product of the particle number and 
the number of the walls of a quasi Neumann condition equals the dependence of the particle 
number of the force correction terms in our model. This implies that the particle number 
dependence can be interpreted as the number of the interactions, as in the case with van der 
Waals equation. 

We also briefly touched upon the three-dimensional box with a quasi-Neumann condition. 
The obtained equation of state (14.191) is a three-dimensional extension of the one-dimensional 
equation (14.161) and we find that the three-dimensional equation possesses the same property 
of the one-dimensional equation, confirming that boundary effects can appear in three di- 
mensions, too. Finally, we mention that generalized boundary conditions have been used for 
effective descriptions of rapidly varying potentials at small scales, which are approximately 
represented by the vanishing limit of widths of step functions [HI dSl ED]- There, the bound- 
ary condition (2.7) and parameter Lq are derived from such limits of two step functions [6]. 
These circumstances may appear, for instance, at junctions of two semiconductor layers and, 
if so, our boundary effects may actually be observed experimentally. 
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A On approximation 

In this appendix, we give a brief account of integral approximation of the sum of series. For 
all n G Z, f(x) can be expanded as 




(A.l) 



where x € (n, n + 1). We obtain 
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where 

oo 1 

T k :=Y,f (k) (n+-). (A.3) 



n=0 



We apply the above equation to f^ 2k \x) and eliminate the differential term by successive 
iteration, then we obtain the following series expansion: 



n qo OO OO 

/ f(x)dx = aste-i/^CO) + f( n + 2- 

^° n=l n=0 



(A.4) 



where 



1 1 



k-i 



d2k-l — 2«2fc+l — 2/C2n+l a 2(fc-n)-l- (A-5) 



n=l 



The coefficient g^-i tends to rapidly in the limit k — > oo. In addition, f(x) = F(y/ax) in 
our system; therefore 

/ (n) (0) ~0[a% (A.6) 

From the above equations, then we can neglect the sum of the first term of the right hand 
side of fTA~4l) if a < 1. 

Now, we show a similar formula to (IA.4D . Let f(x) expand as follows: 

f( x ) = J2v\ x - n ) / (fc) w- ( A - 7 ) 

fc=0 

We define := Y^=of^ k \ n ) an d repeat the above discussion, then the following equation 
is obtained: 

/ - E (^^t/ <6, («) + E E / <2 "(«) 

A/ — — 1 TL — 

OO OO 

= E/( n )+Ek/ (B) (°). ( A - 8 ) 

n=l n=0 

where 

, 1 

°2fc — -«2fc~l 

fc-1 

(A.9) 

n=l 
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B On polylogarithm 



The polylogarithm [22], Li s (z), is defined in \z\ < 1 as follows: 



OO J, 

Z K 



fe=l 

In |z| > 1, it is defined by analytical continuation of the above equation. It satisfies the 
following formulae: 

Li.(-z) + U s {z) = 2 1 - s Li s (^ 2 ) (B.2) 
Jim Li s (e") = -— ^— ; (B.3) 

Re(|i)-+oo 1 [S + 1) 



therefore we obtain the asymptotic behavior of them: 



Lii{ -exp(z)} -> -2y^ (B.4) 
Li f {-ex P (,)}--y^. (B.5) 



We can express Lii(±e y ) and Lis(±e 1 ') in integral forms as 

r°° 2 dr 

Lii(±e») = / — p— (B.6) 

7o ±exp[7re 2 -y\ - 1 

Lis (±e y ) = / = -. B.7 

Jo ±exp[vrx 2 - y] - 1 

The above equations imply that Lii(±e y ) and Li3(±e^) are monotone functions because of 
the integrands are monotone functions of y. 
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